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Abstract. Observational constraints on the primordial deuterium-to- 
hydrogen ratio (D/H) can test theories of the early universe and provide 
constraints on models of big bang nucleosynthesis (BBN). We measure deu- 
terium absorption in high-redshift, metal-poor QSO absorption systems and 
directly infer the value of primordial D/H. We present two measurements of 
D/H, and find D/H = 3.3 ±0.3 x 10~ 5 at z = 3.572 towards QSO 1937- 
1009 and D/H = 4.0 ±0.7 x 10" 5 at z = 2.504 towards QSO 1009±2956. 
Both measurements use multiple-component Voigt profile analysis of high 
resolution, high signal-to-noise spectra and determinations of the Lyman 
continuum optical depth in low resolution spectra to constrain the column 
densities of deuterium and hydrogen. The measurements are consistent with 
a single primordial value of D/H = 3.4±0.3 x 10~ 5 . This is a relatively low 
value, which supports homogeneous models of BBN and standard models 
of galactic chemical evolution. With standard BBN, we find a cosmological 
baryon-to-photon ratio, r\ — 5.1± 0.3 x 10~ 10 , and a present-day baryon den- 
sity in units of the critical density, h\ 00 = 0.019 ± 0.001. These values are 
consistent with high abundance measurements of 4 He in extragalactic H II 
regions|I|, and estimates of 7 Li in warm, metal-poor halo starsJJ. 



1. The Significance of Deuterium 

In the first one thousand seconds of the universe, light nuclei (D, 3 He, 4 He, and 7 Li) 
are created during the epoch of BBN0,|Q|. In the standard modelQ the nuclear yields 
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of the light elements depend on a single parameter, 77, the baryon-to-photon ratio. The 
density of photons in the universe is accurately known from the temperature of the 
cosmic microwave background |J. Therefore, a determination of a single primordial light 
element abundance gives a measure of 77 and the present-day density of baryons, Qb- 
The abundance ratio of deuterium to hydrogen (D/H) is a sensitive function of 77, and a 
measurement of primordial D/H places the strong constraints on 77 and Qb- For example, 
a measurement of D/H with an uncertainty of 20% corresponds to 12% uncertainty in 
both 77 and Qb- 

The determination of the primordial abundance ratio of D/H is a difficult task. 
For the last 25 years, BBN has stood as the only known cosmological source for deu- 
terium nuclei ||. The deuteron is a fragile nucleus, and is easily and totally destroyed 
through stellar processing J?]. With no other source of D besides BBN, any observational 
measurement of D /H will always give a lower limit to the primordial ratio. Local obser- 
vations, in either the solar system or interstellar medium, provide a strong lower bound 
to the primordial abundance, D/H > 1.6 x 10~ 5 ||. Inferring an upper limit from lo- 
cal measurements of D/H requires a problematic extrapolation with models of chemical 
evolution over the last 10 billion years of stellar processing in the galaxy ||. 

In these proceedings, we will describe new measurements of D/H at high redshift 
in intergalactic absorption systems detected along the line of sight to distant quasars 
(QSOs). In 1976, AdamsflCj realized that high redshift systems showing strong Lya 
absorption would provide an ideal site to determine the primordial ratio of D/H. Some 
of the more compelling reasons include (1) the high redshift allows much less time for 
D/H to be significantly altered from the primordial value, (2) the absorption systems are 
very metal poor (less than 1/100 solar metallicity in most cases), which limits the amount 
of destruction of D due to stellar processing, (3) most high redshift absorbers are not 
likely to be associated with galaxies, (4) the systems which are useful in a measurement 
of D/H must have low temperatures (T < 20000 K) and very small turbulent motions, 
which limits the amount of energy which could be introduced into the system. 



2. The QSO Absorption Systems 

Modern spectroscopy of high redshift QSOs reveal hundreds of discrete absorption fea- 
tures which correspond to intervening clouds of gas which lie along the line of sight. A 
small subset of these absorbers can yield a measurement of D/H, the set defined by the 
following criteria: 

The neutral hydrogen column density (number of atoms per square cm) must ex- 
ceed 10 17 cm 2 . 

The velocity structure along the line of sight must be both simple and narrow. 
Simple, in that only a few components with small velocity separations are detected 
in absorption. 

The redshift of the absorber must be high, z a b s > 2.5, to shift all of the desired 
absorption features into optical wavelengths which can be observed by large ground 
based telescopes. 
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In addition, we prefer systems which show weak or undetectable metal absorption lines, 
and lie along the line of sight to the brightest QSOs (which makes high-resolution spec- 
troscopic observations practical). 

The criteria listed above places severe limits on the sample of absorption systems. 
Out of the hundreds of absorbers along the line of sight to each QSO, we find systems 
to measure D/H in less then 3% of the lines of sight. That is, less than 0.01% of the 
absorbers in the entire observable sample are suitable systems. Needless to say, QSO 
absorption systems which show deuterium are rare. 

The recent success of the high redshift D/H measurements rely on a method to 
increase the likelihood of finding these rare systems. Systems with the required high 
H I column densities are optically thick to photons with energies above one Rydberg, 
and exhibit very prominent continuous absorption in QSO spectra. These systems are 
known as Lyman limit systems. Lyman limit systems are easy to detect in low resolution 
spectra, which gives a much faster means to find likely systems, compared to obtaining 
high resolution data of random lines of sight to look for deuterium absorption. The top 
panel of Figure 1 shows a low resolution spectrum of QSO 1937-1009. The Lyman limit 
system at z = 3.572 is responsible for the continuous absorption of flux below 4200 A. 
The Lyman limit serves not only to identify the systems, but the optical depth of the 
absorption is directly proportional to the H I column density in the system. The role of 
the Lyman limit measurements will be discussed in detail in the next section. 

Figure 1 shows a schematic view of two absorption features in the deuterium system 
towards QSO 1937-1009. The low resolution spectrum in the top panel covers a large 
spectral range and is plotted in units of flux per A. The spectrum shows the intrinsic 
shape of the QSO continuum, including emission features, but also depicts the sharp 
increase in absorption at wavelengths shorter than Lya emission at 5850 A, due to the 
intervening Lya absorption systems (the "Lya forest"). The next set of panels zooms 
in on two specific regions in the spectrum, corresponding to the wavelengths of Lya and 
singly ionized carbon (C II) redshifted to z — 3.572 (A = Ao(l + z)). The two regions 
are noticeably different. The left hand panel shows the high density of lines in the Lyct 
forest at high redshift surrounding the strong Lya line of the deuterium system. The 
right hand panel shows only the weak optically thin absorption from C II, and shows 
the absence of any other features nearby, so there is no confusion. The bottom panels 
show the profiles of the Lya and C II. Deuterium absorption can be see on the left wing 
of the H I Lya feature. The isotopic separation of D and H lines is -82 km s _1 ; the D I 
energy levels are shifted up in energy by a factor of 1.000272 compared to hydrogen, 
due to the small change in the reduced mass from the additional neutron. The middle 
and bottom panels show spectra obtained with the high resolution echelle spectrograph 
(HIRES) |n]]. The spectra have a resolution of 8 km s _1 FWHM, which is sufficient to 
resolve the deuterium lines. The spectra have been normalized to allow a direct analysis 
of the absorption features. The intrinsic QSO continuum was modeled with a low order 
polynomial in each region. The next section will describe the method to measure D/H 
from the Lyman series absorption profiles of D I and H I. 

3. The Deuterium Measurements 

Here we present the spectra and analysis of the two absorption systems which yield 
the best measurements of D/H at high redshift. For a more detailed description of the 



3 



01 



o 20 



g? 10 

CO 

2 o 

x 



0.5 



x 

£ 

CD 

N 
• i— i 

a i 

o 



~i 1 1 r 



n 1 1 r 



~i 1 1 r 



~i r 



QSO 1937-1009 



J abs 



3.572 




5540 / 556Q\5580 6080 /6100X6120 



0.5 - 







D H 






C II (1334)- 

I I I I I I I I I I I I I I I I ~ 



5555 



5560 6101 
Wavelength (A) 



6103 



Figure 1. QSO 1937-1009: See text for more details 
Top: Lick spectrum (FWHM = 4 A) 

Middle and Bottom: Keck+HIRES spectrum of the Lya region (left) and C II region 
(right) at z = 3.5722 
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method and analysis, the reader is referred to larger papers[12, 13, 14, 15, 16 
description is given here. 



A brief 



3.1. The Method to Measure D/H 

We construct a model fit with discrete absorption lines, given by Voigt profiles, convolved 
with the instrumental resolution. Each Voigt profile is given by three parameters, the 
column density (N), the redshift (z), and the intrinsic velocity dispersion along the line of 
sight (b). Each hydrogen line profile with a large hydrogen column, N(H I) > 10 16 cm -2 
has a corresponding deuterium profile at the same redshift. The velocity dispersion of 
both the H I and D I profile are given by two parameters, the temperature and turbulent 
velocity dispersion. In short, we substitute the free parameters of b(H I) and b(D I) 
with two parameters, T and b tur , which are identical in both H I and D I absorption 
lines. We also assume that D/H is global, that is, all components have the same ratio of 
N(DI)/N(HI). 

We select regions of the HIRES spectrum which contain the absorption features of 
the Lyman series and are not severely blended by unrelated lines in the Lya forest. The 
continuum is modeled as a low order polynomial in each region, and the coefficients of 
the polynomial are treated as free parameters in the fitting process to statistically allow 
for continuum uncertainties. 

In this analysis, a measure of D/H is the goal, so we study the goodness of fit x 2 
as a function of D/H. To directly study the effect of D/H on the model fit, we model 
the spectra with hundreds of iterations. In each iteration, we specify D/H, let all other 
parameters go free, and find the best x 2 associated with the specified value of D/H. 
We iterate over a large range in D/H, and measure the most likely value of D/H by 
locating the minimum of x 2 as a function of D/H. This method has much to offer. The 
parameters describing the profiles and continuum are sometimes highly correlated. The 
true uncertainty in D /H can not be determined by calculating the individual parameter 
uncertainties near the minimum of the x 2 function. We must explicitly map the x 2 
dependence on D/H to include the correlations which are present in the model. 



3.2. QSO 1937-1009 

Tytler et al (1996) |]15] made the first measurement of low D/H in the absorption system 
at z = 3.572 towards Q1937-1009. We analyzed the high-resolution spectrum (8 hrs 
of exposure), which resolved the entire Lyman series up to Ly-19, as well as associated 
metal lines. By profile fitting the Lyman lines, with the position of the velocity com- 
ponents given by the metal lines, we find D/H = 2.3 ±0.3 ± 0.3 x 10~ 5 (statistical and 
systematic errors). The largest uncertainty in the measurement is the neutral hydrogen 
column density, log N(H I) = 17.94 ±0.06 ± 0.05, and the uncertainty stems from the 
saturated Lyman profiles (discussed in detail below). We then obtained a high quality 
low-resolution spectra from Keck with LRIS|T7|], which gave better sensitivity short- ward 
of the Lyman limit, to directly measure the total N(H I) in the system and therefore 
place better constraints on D/H. Utilizing both the high and low-resolution spectra, we 
find log N(H I) = 17.86 ±0.02 by a direct measurement of the optical depth short-ward 
of the Lyman limit at 4200 A (Buries & Tytler 1997a). With this constraint and the 
more sophisticated method described above, we measure D/H = 3.3 ± 0.3 x 10 -5 (Buries 
& Tytler 1997b). 
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Figure 2. Lyman series of Q1937-1009 at z = 3.5722. 
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Figure 3. Lyman limit region of Q1009+2956 at z — 2.504. 

Top: HIRES spectrum (FWHM = 8 km/s) showing lines Ly-11 to Ly-24. Bottom: Lick 
spectrum (FWHM = 4 A) and model fit short-ward of Lyman limit. The solid line shows 
the best fit and dotted lines show la errors. 
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Figure 4. Lya region of z — 2.504 system towards Q1009+2956. The histogram shows 
the pixels of the HIRES spectrum normalized to the initial continuum estimate (solid 
line at unity). The model profile is composed of the main H I (dot-dashed) and D I 
(dashed), and contaminating H I (dotted). 



The Lyman lines used in the fitting procedure are shown in Figure 2. The spectra 
have been shifted to a velocity scale by the simple transformation, - c = Xo ^ +z ^ — 1-0, 
where Ao is the rest wavelength of the Lyman line, and z is the redshift corresponding 
to zero velocity, in this case z = 3.57221. 

Although the positions of the metal line components were previously used to con- 
strain the fit 1 15], the present analysis does not require any information from the metal 
lines in the fitting procedure. The use of the metal lines can be justified for this system, 
but in general, one is likely to introduce systematics by including the metal lines, and it 
remains a good policy to leave all other lines except H and D out of the fitting procedure. 



4. QSO 1009+2946 

We discovered an absorption system at z — 2.504 towards Q1009+2956 ideal for a 
measurement of D/H(16[. This system has a lower hydrogen column density, log N(H I) 
= 17.39 ± 0.06. The highest order Lyman lines become unsaturated, which yields a 
precise measurement of N(H I) in both low and high resolution spectra (Fig. 3). Over 
twelve hours of Keck+HIRES produced a very high quality spectrum of the entire Lyman 
series, resolving the entire series up to Ly-22. We find strong evidence for contamination 
of the deuterium Lya absorption feature, which introduces the largest uncertainty in 
the measurement. Figure 4 shows the Lya line responsible for the contamination. The 
profiles of D I (dashed lines) and H I (dash-dotted lines) are shown explicitly, with the 
weak contaminating hydrogen line (dotted) falling just blueward (to the left) of D I. We 
include this contamination as a free parameter and find D/H = 4.0 ±0.7 x 10 5 1|1 3|| - 
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5. Conclusions 



These two measurements of D/H in QSO absorption systems are the best and most 
robust measures to date. Deuterium has been identified and analyzed in a number of 
other QSO absorption systems [[18f] We have found another two systems which place a 
strong upper limit on D/H at D/H < 10~ 4 . Combined with the two measurements 
described above, the four independent systems support a low primordial abundance of 
deuterium, and together give D/H = 3.4 ±0.3 x 10 -5 . If this represents the primordial 
value, nucleosynthesis calculations from standard BBN models with three light neutrinos 
give 77 = 5.1 ± 0.3 x 10" 10 and tt b h\ m = 0.019 ± 0.001. 

The constraints from D/H can be utilized to constrain cosmological models, quan- 
tify dark matter both in unobserved baryons and non-baryons, specify the zero point 
for models of deuterium evolution]^], test directly the predictions of standard BBN by 
comparing with other light element abundances |19[, and limit the amount of small scale 
entropy fluctuations in the early universe [I 
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